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ABSTRACT 

The binding energy parameter A plays an important role in common enve- 
lope (CE) evolution. Previous works have already pointed out that A varies 
throughout the stellar evolution, though it has been adopted as a constant in 
most of the population synthesis calculations. We have systematically calculated 
the binding energy parameter A for both Population I and Population II stars of 
masses 1 — 20Mq, taking into account the contribution from the internal energy 
of stellar matter. We present fitting formulae for A that can be incorporated into 
future population synthesis investigations. We also briefly discuss the possible 
applications of the results in binary evolutions. 

Subject headings: stars: evolution - stars: mass loss - binaries: general 
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Introduction 



One of the key stages in the evolution of close binary stars is the common envelope 
(CE) evolution. When the primary star fills its Roche-lobe (RL) and the following mass 
transfer is dynamically unstable so that the secondary cannot accrete all the transferred 
material, the transferred matter will form an envelope embedding both stars and the 
binary enters the CE evolution. The secondary star orbits inside the CE, and experiences 
a drag force by the envelope. This causes the or bital decay and spiral-in of the star wit h 



possi ble ejection of the envelope (for reviews see 



Iben &: Livio 



1993 



Taam &: Sandquist 



20001 ). Two main kinds of theories have been developed to describe the envelope ejection 



process. One is the a formalism, whic h considers the conversion o 



overcome the envelope binding energy (jWebbink 



1984 



de Kool 



1990 



orbital energy to 



Dewi fc Tauris 



20001 ). 



If the orbital energy of the secondary is large enough to eject the envelope, the system will 
survive and become a compact binary containing the core of the lobe-filling star and the 
secondary in a much smaller orbit; if the orbital energy is not enough, the two stars will 
coalesce. The other is t he 7 formalism, which co nsiders angular momentum transformation 



during the spiral-in (see 



Nelemans fc Tout 



2005 



and reference therein). Despite extensive 



three dimensional 



lydrodynamical simulations of CE evolution (for recent reviews, see 



Taam fc Sandquist 



2000 



Taam fc Ricker 



20061 ). the physics of the CE phase remains 



poorly understood, including the efficiency with which the CE is ejected from the system. 
Al ternative evolutionary pat hs h ave been proposed to avoid the CE evolution, such as those 



by 



King fc Begelmanl (119991 ) and 



Beer et al. 



(I2OO7I ). who suggested that super-Eddington 



accretion may lead to mass ejection from the system instead of entering the CE stage. 



According to iHan et al. 



(11994 



19951) and IPewi fc TaurisI (I2000h . the total binding 



energy of the envelope is described by: 



bind 



dm + ath 



Mc, 



Udm, 



(1) 
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where Mdonor is the mass of the donor star and Mcore is its core mass, —GM{r)/r and U are 
the gravitational and internal energy of the stellar matter respectively, G is the gravitational 
constant, and ath the percentage of the internal energy contributing to ejection of the 
envelope, normally taking the value between and 1. 



For conven i ence the tota l binding energy of t 



flWebbink 



1984 : 



de Kool 



1990 



Dewi &: Tauris 



le envelope is usually expressed as 



2000|): 



Si 



C Mdonor ^env 



bind 



AoirL 



(2) 



where M^nv is the envelope mass, A is the binding energy parameter, tl = Rl/O'I is the ratio 
of the RL radius and the orbital separation at the onset of CE, and oirL is normally taken 
as the stellar radius once a star fills its RL and starts to transfer matter. If we assume 



that part of the change of the orbital energy is used to eject the envelope, as described by 
3 



WebbinkI (ll98J), then 



Eund — OlCE 



(3) 



2af 2ai 

where acE is the efficiency of converting the orbital energy to kinetic energy to eject the CE, 
M2 the mass of the secondary, and ai and Of refer to the initial and final orbital separation 
of the CE phase, respectively. 

Combine Eqs. (2) and (3), the final orbital separation is related to the pre-CE 
separation by the following formula. 

Of ^ Mco,eM2 1 

fli Mdonor M2 + 2Menv/acEArL' 

Equation (4) indicates that the orbital separation after CE is sensitively dependent on the 
product of acE and A, and the convenient way is to treat both acE and A as constants of the 
order of unity. In most of the population synthesis investigations concerning the evolution 
of close binaries in the literature, A was usually taken to be a constant, normally around 0.5 
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throughout the evolution o : 



1994 



Dewi fc Tauris 



A is a variable. 



200d. 



the system. However, there have lots of works (e.g. iHan et al 



2001 



Dewi k TaurisI fl2000h and 



Podsiadlowski et al. 



2003 



Webbink 



20071 ) suggesting that 



Dewi &: TaurisI ( 1200 ll ) found that the value of A 



changes as the star evolves, and reaches far more than 0.5 in late evolutionary stages for 
stars with mass b etween 3 and QMq. In recent population synthe sis on post-CE b inarie s, 



Davis et al. 



( I2OIOI ) used line ar interpolation s of th e tab ular data in iDewi fc 



calculate A. We notice that 



Dewi fc TaurisI mm and 



Dewi &: TaurisI (l200ll )'s calculations 



aurisi mm to 



did not cover Pop. 1 stars less massive than ?)Mq and Pop. 11 stars, and their results 
in tabular form may not be easily used. Obviously more advanced population synthesis 
requires a simple description of the parameter A for better constraints on the initial stellar 
population. In this work we are attempting to investigate the binding energy parameter 
A by considering the contribution of internal energy of stellar matter for both Pop. I and 
Pop. II stars, and to find easy-to-use fitting formulae of A for further research. 

This paper is organized as follows. We describe the stellar evolution code and 
assumptions adopted in Sect. O In Sect. [3] we present the calculated results. We summarize 
our results and briefly discuss their possible imphcations in Sect. jH 



2. Model Description 



We adopt th e latest st e. 
modified later by 



Egdetoni (Il972l l 



lar evolution code TW! 



Eggleton et al. 



mm and 



develo ped bv lEggletonl (Il97ll) and 



Pols et al. 



(119951 ). One of 



^The code used here is a variant of the code ev described, in its initial version, by Eggleton 
(1971, 1972) and Eggleton et al. (1973), and was developed by the present authors from the 
version of Eggleton (2009, unpublished). The current version of ev is obtainable on request 
from eggletonl@llnl.gov, along with data tables and a user manual. 
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the major improvements of the code is that, when a low-mass star is reaching the He flash 
igniting point, the code will automatically generate a zero- age horizontal branch star with 
both the same core and total mass which enables further evolution. 

We take the mixing length parameter a = I /Hp to be 2.0 and follow the evolution 
of both Pop. I and Pop. II stars. We assume the chemical composition of X = 0.7 and 
Z = 0.02 for Pop. I stars and X = 0.755 and Z = 0.001 for Pop. II ones. Once a star 
evolves off the main- sequence, we use the layer that contains less than 15% H as the 
boundary between the core and the envelope. We then calculate the envelope binding 
energy from Eq. ([T]) by setting ath to be either or 1. After that we use Eq. ([2]) to 
calculate the corresponding value of A, where we simply take aiVi^ = Ri^ as the stellar radius 



R. To include the stel 
mechanisms similar to 



ar wind mas s loss during the evolution, we consider the following 



Hurley et al. 



(I2OOOI). 



1. The modified Nieuwenhuijzen-type wind model (INieuwenhuijzen &: de Jageiill990l ): 

Z 



NJ 



9.6 X 10-i5(— )i/2^°-*^i^i.24^o.i6^^y^- 

Zq 



(5) 



where M, i?, and L are the stellar mass, radius and luminosity in solar units, respectively. 

OSS r ate is modified by the factor Z^/'^ due to different metallicity 



The wind mass 


OSS r 


(Kudritzki et al. 


1987) 



from the fact that the effective number of lines contributing to the line acceleration changes 
with metallicit y Z and, he r ice, c auses a metallicity dependence of the mass-loss rate. The 



observations of 



Puis et al. 



(Il996[ ). which in the metallicity range between the Galaxy and 



the SMC {Z ~ Zq/^) confirms the exponent of 1/2 for 0-stars. 



2. The Reimers type wind fo r giant-branch stars as described by 



Kudritzki fc Reimers 



f!1978h and llben k Renzinil fll983l ): 



Mj, = 4 X 10-13^ 
M 



(6) 
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with T] = 0.5 to modify the wind loss. 



and 



3. The reduced Wolf-Rayet-l ike mass loss by iHamann. Koesterke fc Wessolowskil (Il995[ ) 



Hamann &: Koesterke 



([19981): 

MwR = 10-'^L' \1.0 - fi)MQjT-\ (7) 

with 

/i= (^^^^)min{5.0,max[1.2,(A)^]J (g) 

where Lq = 7.0 x 10'^ and k = —0.5. 

I 



4. AGB stellar wind loss bv iVassiliadis fc WoodI (11993[ l: 



logMvw = -11-4 + 0.0125[Po - 100max(M - 2.5, O.O)]M0yr-\ (9) 

where 

log Po = min(3.3, -2.07 - 0.9 log M + 1.94 log R), (10) 
and this wind loss is hmited to 

Mvw,max = 1.36 X lO-^LMeyr-^ (11) 

Then the total mass loss via stellar wind (absolute value) is taken to be 

M = max(MNj, Mr, Mwr, Mvw)- (12) 

When taking into account the contribution of internal energy to total binding energy, 
we include the following terms: the thermal energy of the stellar matter, the radiation 
energy, the ionization of H, H"*", He, He"*" and He"'""'", the dissociation of H2, the ionization of 
C, N, O, Ne, Mg, Si and Fe (these seven elements are assumed to be fully ionized for all the 
temperatures). To descriminate the effect of internal energy, we set ath to be and 1, and 
calculate the respective binding energy parameter, namely Ag and Ab, respectively. 
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3. Calculation Results 

We have made calculations for both Pop. I and Pop. II stars, with mass range from 

1 to 20 Mq. For each star we follow its evolution until its age exceeds 15 Gyr or the 
code crushes for rapid C burning. Once a star evolves to produce a H-exhausted core we 
calculate the binding energy of its envelope and accordingly Ab and Ag, both with and 
without the internal energy of the stellar matter. In Figs. [1] and |2] we show some examples 
of the evolution of Ab and Ag with respect to the stellar radius R. 

Since a star experiences both expansion and shrinkage throughout its life, we divide its 
evolution into three stages after the star leaves the main-sequence according to the change 
of the stellar radius, to make the fitting of A more practicable. Stage 1 begins at the center 
H exhaustion and ends when the star starts to shrink (i.e., near center He ignition). Stage 

2 follows and ends when the star starts to expand agaiij^. Stage 3 begins after that and 
continues till the end of the evolution. 

Figures [1] and [2] show that both Ab and Ag vary during the evolution of stars. For IMq 
star, both Ab and Ag decrease with i?, but their magnitudes are around unity before the 
He flash. After that the values of A have a big decrease along with the stellar radius, but 
Ab ~ 2Ag throughout the evolution. More massive stars often experience rapid increase 
in A in stage 3. For stars with mass from ~ 3 to SM© (the mass range is related to 

^Generally a star shrinks after its center He ignition and starts to expand again later. 
For convenience of fitting we define the end point of stage 2 (and also the starting point 
of stage 3) as the minimum of stellar radius. For example, for Pop I stars which are more 
massive than lOM©, we use 10^'^ center He abundance as the end point of stage 2. As a less 
massive star expands again before its center He exhaustion, we simply assume that stage 2 
ends when the respective HB star starts to expand again. 
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metallicity), Ab and Ag take the values ~ 1 before the star reaches stage 3, during which 
the internal energy dominates, and the total binding energy finally becomes positive. The 
resulting Ab increases rapidly and eventually turns out to be infinity, while Ag remains to 
be around 0.5. For stars massive than 6M0, the binding energy decreases in stage 3 but 
never becomes positive, correspondingly Ab takes smaller value and Ag is around 0.5. For 
stars more massive than 12Mq, A even drops below 0.5 for most R. The extreme case in 
our calculation is for 20Mq stars, in which both Ab and Ag evolve to be -C 1 in stage 3. 

To provide easy-to-use formulae of Ab and Ag, we make polynomial fitting with the 
following formulation, 

y = a + bix + h2X^ + h^x'^ + h^x'^ + fesx^, (13) 

where y is Ab or Ag, and x is the stellar radius R in most cases. For low- mass stars, we 
also use the fractional r nass o f the envelope m^^^ = Menv/Mdonor as the fitting variable as 



suggested by IWebbinkI (j2007l ). In limited cases we use logAb instead of Ab as the y variable 
during stage 3 when the value of Ab is too large, and take 1/X as y and mgnv as x in several 
cases (see Tables [1] and [2] for details). In Tables [1] and [2] we present the fitting parameters 
for Pop. I and Pop. II stars, respectively. 



4. Discussion and conclusions 

Our calculations show that stars with different masses have different values of A, and 
A is not constant for the same star in different evolutionary stages. From Figs. [H and [21 Ab 
diverges from 0.5 for most M and R, and it is obvious that assuming a constant A = 0.5 
in population synthesis calculations is far from fact and therefore lack of reliability. It is 
also interesting to note that the range of Ag is much narrower than that of Ab. Of course 
the actual value of A should lie between Ag and Ab since not all of the internal energy 
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contributes to the ejection of the envelope (jPewi fc Taurial2000l ). As we have given the 
fitting formulae for A, the results can be useful for future population synthesis works. 



However, this approach shou 
binding energies are positivqf 



d be adopted carefully, especially in the cases when the 
and where the efficiency of conversion of ionization energy 
to kinetic energy may differ from the efficiency of conversion of other forms of internal 
energy. As we can see from Eqs. ([T]) and the calculated values of A depend on the 
stellar mass, the core mass, the envelope mass distribution and the percentage of internal 
energy that contributes to envelope ejection. Factors like metallicity and stellar wind mass 
loss also affect the stellar evolution and the resulting values of A. Additionally, the different 
definitions of core-envelope boundary lead to different stellar core mass and therefore affect 
the calculated results. 



and 



Our results agree rea sonably with that of 



HanetaL 



fll994l ). bewi fc TaurisI fcoooh 



Dewi &: TaurisI ( 120011 ). although there exist some differences. 



Dewi fc TaurisI (120001 ) 



proposed that A is almost independent of the stellar chemical composition. From our 
calculations we find that given the same mass. Pop. I and Pop. II stars have different 
A-values (see Figs. [T]and[2]). For example, the mass range in which the values of Ab can 
become negative is narrower for Pop. II stars than for Pop. I ones. In Fig. [3] we compare 
the evolution of the two As for a 5Mq star with different metallicity. It is seen that the A 
evolutions show little difference at the beginning, but later they diverge from each other, 
and the values of A for the Pop. I star always lie above those of the Pop. II star in late 



■^Although the gravitational binding energy remains negative for all the stars throughout 
their evolution, the total binding energy, on the other hand, can turn to be positive for stars 
with mass ranging from ~ 3 — 5Mq (specifically, 3 to 5.6Mq for Pop. I stars and 2.8 to 
4.5M0 for Pop. II stars at the end of their AGB phase when their radii expand to be more 
than hundreds of solar radius, and ionizaiton energy dominates the total binding energy). 
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evolutionary stages. The reason is as follows. The binding energy (absolute value) of the 
Pop. I star is either comparable with (in early evolutionary stage) or lower than (due to 
the contribution of ionization energy) that of the Pop. II star. Since the Pop. I star evolves 
slower than the Pop. II one, as a result the Pop. I star has a smaller H-exhausted core and 
higher-mass envelope than the Pop. II star when they have the same radii, thus higher A 
values. 

We used the 15% H abundance layer as the stellar core-envelope boundary to calculate 
the binding energy and A, since the core mass from such definition increases in most of 
the time and is suitable for calculating the binding energy (Eggleton 2009). Note that 



Dewi fc TaurisI ( 120001 ) used the 10% H abundance layer as the core-envelope boundary, and 
in our calculation we found that the two definitions give similar results. There are other 
criteria to define this boundary, such as the maximum energy generati on layer or the density 



gra dient criteria (iHan et al. 



and 



1994 : 



Dewi &: Tauris 



20001 ). According to 



Tauris fc Dewil fl200l[ ) 



Dewi &: TaurisI (l200l[ ). most of these criteria give outer core boundary and higher core 



mass. In general, if the total stellar mass is a constant, a more massi ve core will decrease 



the total binding energy of the enve 



Tauris fc Dewi 



2001 



ope a nd result in a larger A- value ( iDewi fc Tauris 



Dewi fc Tauris 



2000 



20011 ). In other words, the A- values we calculated may 



be regarded as the lower limit of the actual ones. Stars more massive than IOMq lose 10% 



to 30% of it s total mass during its evo 



valu es of A ( 
and 



Dewi &: Tauris 



Podsiadlowski et al. 



2000 



2001 



ution, which can affect the s tellar structure a.nd th e 



Podsiadlowski et al. 



20031 ) 



Dewi fc TaurisI (1200 If ) 



(120031 ) have also found that A is sensitive to the wind mass loss for 
massive stars, and including stellar wind leads to higher binding energy (absolute value) 
and lower value of A. 

Our results suggest that Ab can take values significantly larger than 0.5 during the 
late-stage evolution of stars with mass between ~ 2 and IOMq, which will significantly 
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increase the efficiency of envelope ejection. If we consider an evolved giant or AGB star in 
this mass range with A >> 0.5, the envelope binding energy may be much smaller compared 
to that with A = 0.5. It is then much easier for its companion to eject the CE, and to 
produce wider binaries after the CE phase. It is also worth noticing that both Ag and Ab 
are less than 1 in the late evolutionary stages of massive stars. For stars massive than 
~ 8Mq, a even lies below 0.2. As a result, the actual envelope ejection efficiency for these 
stars decreases considerably. Taking a 2OM0 Pop. I giant star as an example. As seen from 
Fig. [H both Ab and Ag are <^ 1, which means that the envelope binding energy become 
much higher. The secondary has to convert more orbital energy to eject the envelope, which 
probably leads to a much tighter final orbit, or even a merger. These may have important 
implications for the formation and evolution of related binaries that produce compact stars, 
which will be discussed in future works in more detail. 

We thank Prof. P. P. Eggleton for providing the latest version of the evolutionary code 
and helpful suggestions. This work was supported by the Natural Science Foundation of 
China (under grant number 10873008) and the National Basic Research Program of China 
(973 Program 2009CB824800). 
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Fig. 1. — The evolution of the binding energy parameter A for Pop. I stars of mass IMq to 
20Mq. The upper and lower solid lines are Ab and Ag, respectively. The dot-dashed lines are 
the fitting results. 
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Fig. 2. — Same as Fig. 1 but for Pop. II stars. 
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Fig. 3. — Comparison of the A values for SMq Pop. I and Pop. II stars. The upper and 
lower solid lines represent Ab and Ag for the Pop. I star, respectively. The upper and lower 
dotted lines are Ab and Ag for the Pop. II star, respectively. 



Table 1. Fitting Results of A for Pop. I Stars 



Stage Mass (Mq) A a 61 62 63 64 65 



1** 




19.84809 


-69.00242 


81.24502 


-31.5572 








** 


As 


26.09283 


-72.37166 


64.13808 


-16.24247 








2 


Ai, 


2.05363 


-0.00685 


-3.42739E-4 


3.93987E-6 


-1.18237E-8 







As 


1.07658 


-0.01041 


-4.90553E-5 


1.13528E-6 


-3.91609E-9 





3 


Ab 


2.40831 


-0.42459 


0.03431 


-9.26879E-4 


8.24522E-6 







As 


1.30705 


-0.22924 


0.01847 


-5.06216E-4 


4.57098E-6 





4 


Ai, 


1.8186 


-0.17464 


0.00828 


-1.31727E-4 


7.08329E-7 







As 


1.02183 


-0.1024 


0.00493 


-8.16343E-5 


4.55426E-7 





5 


A6 


1.52581 


-0.08125 


0.00219 


-2.0527E-5 


6.79169E-8 







As 


0.85723 


-0.04922 


0.00137 


-1.36163E-5 


4.68683E-8 





6 


A(, 


1.41601 


-0.04965 


8.51527E-4 


-5.54384E-6 


1.32336E-8 







As 


0.78428 


-0.02959 


5.2013E-4 


-3.45172E-6 


8.17248E-9 





7 


At 


1.38344 


-0.04093 


5.78952E-4 


-3.19227E-6 


6.40902E-9 







As 


0.76009 


-0.02412 


3.47104E-4 


-1.92347E-6 


3.79609E-9 





Q 

cS 


A(, 


i.OOOiu 






1 QC^QQI TT K 

-i.aoyoiiij-D 


o.yjoooziij-tf 


u 




As 


0.73826 


-0.01995 


2.37842E-4 


-1.09803E-6 


1.79044E-9 





9 


A(, 


1.32549 


-0.02845 


2.79097E-4 


-1.07254E-6 


1.46801E-9 







As 


0.71571 


-0.01657 


1.64607E-4 


-6.31935E-7 


8.52082E-10 





10 


A(, 


1.29312 


-0.02371 


1.93764E-4 


-6.19576E-7 


7.04227E-10 







As 


0.69245 


-0.01398 


1.17256E-4 


-3.81487E-7 


4.35818E-10 





12 


A6 


0.23031 


0.74955 


13.20774 


0.63332 


68.67896 







As 


0.10247 


0.31709 


60.57985 


0.7254 


8.90151 





14 


Ab 


0.0576 


0.59416 


158.80613 


0.87013 


19.04369 







As 


0.07178 


0.29286 


80.4363 


0.73414 


10.52307 





16 


Ab 


-0.2064 


0.88945 


22.23505 


0.79813 


287.58284 







As 


0.04116 


0.74011 


11.60388 


0.29533 


97.68433 





20 


Ab 


0.04959 


0.85626 


29.90918 


0.44041 


199.20324 







As 


0.04101 


0.24106 


105.06077 


0.68656 


14.81172 







A(, 


16.71634 


-26.38217 


-14.18029 











** 


As 


24.27819 


-31.47029 


-33.30049 
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Table 1 — Continued 



;age 


Mass (Mq) 


A 


a 


6i 


62 


63 


64 


6b 


2 


2 


A6 


34.41826 


-6.65259 


0.43823 


-0.00953 














13.66058 


-2.48031 


0.15275 


-0.00303 








2 


3 


A6 


-42.98513 


7.90134 


-0.54646 


0.01863 


-3.13101E-4 


2.07468E-6 






Ag 


-6.73842 


1.06656 


-0.05344 


0.00116 


-9.34446E-6 





2 


4 


A6 


-7.3098 


0.56647 


-0.01176 


7.90112E-5 












Ag 


-3.80455 


0.29308 


-0.00603 


4.00471E-5 








2 


5 


A6 


-9.93647 


0.42831 


-0.00544 


2.25848E-5 












Ag 


-5.33279 


0.22728 


-0.00285 


1.16408E-5 








2 


6 


A6 


13.91465 


-0.55579 


0.00809 


-4.94872E-5 


1.08899E-7 









Ag 


7.68768 


-0.30723 


0.00445 


-2.70449E-5 


5.89712E-8 





2 


7 


A6 


4.12387 


-0.12979 


0.00153 


-7.43227E-6 


1.29418E-8 









As 


2.18952 


-0.06892 


8.00936E-4 


-3.78092E-6 


6.3482E-9 





2 


8 


A6 


-3.89189 


0.19378 


-0.0032 


2.39504E-5 


-8.28959E-8 


1.07843E-10 






Ag 


-2.24354 


0.10918 


-0.00179 


1.33244E-5 


-4.57829E-8 


5.90313E-11 


2 


9 


A6 


0.86369 


-0.00995 


4.80837E-5 


-6.10454E-8 


-2.79504E-12 









Ag 


-0.7299 


0.0391 


-5.78132E-4 


3.7072E-6 


-1.07036E-8 


1.14833E-11 


2 


10 


A6 


0.74233 


-0.00623 


2.041971ij-5 


-1.30388Ji-8 












As 


0.36742 


-0.00344 


1.27838E-5 


-1.0722E-8 








2 


12** 


Ab 


-7261.84589 


43423.48906 


-96937.03708 


95946.58045 


-35560.07607 







** 


As 


-280.84015 


1259.61941 


-1674.61224 


681.20633 








2 


14** 


A6 


-3334.04459 


20638.88263 


-47238.6357 


47684.78906 


-17966.67435 







** 


As 


-4769.9922 


29195.29936 


-65789.75807 


65234.39892 


-24107.00589 





2 


16** 


A6 


-321.34016 


1789.56304 


-3002.4542 


1600.03159 










** 


As 


-683.28449 


3792.87458 


-6376.86784 


3419.38525 








2 


20** 


Ai, 


72.31855 


482.85039 


-2593.27601 


3880.69092 


-1913.27392 







** 


As 


166.95028 


649.50464 


-3805.22907 


5447.28964 


-2487.34665 





3 




Ab 


75.23035 


-763.8397 


3180.89498 


-4665.23406 










** 


As 


121.25389 


-1169.84175 


4686.7113 


-6749.8215 








3 


2 


A6 


0.88954 


0.0098 


-3.1411E-5 


7.66979E-8 












As 


0.48271 


0.00584 


-6.22051E-5 


2.41531E-7 


-3.1872E-10 
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Table 1 — Continued 



Stage Mass (Mq) A a 61 62 63 64 65 



Q 


Q* 






U.UU 1 










n m 1 09 

U.Ul iUZ 


Q 


A* 




A Q7Q99 








\ 


n 1 Q1 CQ 
U. lOJ-OO 


u.uuyc5'4: 


Q 


c* 




n anni i 


n nnQQ9 
u.uuyyz 






\ 

'^S 


-U.UU^OD 




Q 









U.UD4D / 






^g 




-U.UUZDD 


3 


7 


Ab 


-0.63266 


0.02054 








0.26294 


-0.00253 


3 


8 


A6 


-0.1288 


0.0099 






Ag 


0.26956 


-0.00219 


3 


9 


A6 


1.19804 


-0.01961 






Ag 


0.40587 


-0.0051 


3 


10 


A6 


0.3707 


2.67221E-4 






Ag 


0.25549 


-0.00152 


3 


12 


A6 


-58.03732 


0.23633 






Ag 


-15.11672 


0.06331 


3 


14 


A6 


-106.90553 


0.36469 






Ag 


-39.93089 


0.13667 


3 


16 


A6 


-154.70559 


0.46718 






Ag 


-65.39602 


0.19763 


3 


20** 


A6 


-260484.85724 


4.26759E6 




** 


Ag 


-480055.67991 


7.87484E6 



A Q979fi1T ^ 
-^L.OZ t ZOII/-O 


Q1 0/1 917 Q 


U 


U 


« /1(^R90TT' ^ 
-D.4D0zyr!j-0 


O.DDoO / rj-y 


7 91 Q1 QTT 1 n 


1 99ni 1? 1 9 


-0 . ZOUOO-Cj- 


K Q7Q9Q'P Q 
0.0 ( oZO-Cj-O 


u 


n 
U 


9 aQaQ9Tr 
- z . csy cSoZxIj- 


9 KQC^I QTT Q 

z.Dooiyxij-o 


u 


U 


Q nQ9/i7"l7 c 
-o.UoZ^ i Sli-O 


O.ZOZ0OI1/-0 


u 


n 
U 


/I 71 1 1 71? K 
4t. i ill i Jzj-D 




u 


U 






U 


U 


9 91 7QQT? 


-Z.ooDyDilj-o 


U 


U 


-1.3646E-4 


2.8661E-7 








1.32272E-5 


-7.12205E-9 








-6.71455E-5 


1.33568E-7 








7.97743E-6 


-1.53296E-9 








1.28222E-4 


-3.41278E-7 


3.35614E-10 





2.73866E-5 


-5.74476E-8 


4.90218E-11 





-9.86464E-6 


2.26185E-8 








3.35239E-6 


2.24224E-10 








-3.20535E-4 


1.45129E-7 








-8.81542E-5 


4.0982E-8 








-4.1472E-4 


1.57349E-7 








-1.55958E-4 


5.94076E-8 








-4.70169E-4 


1.57773E-7 








-1.99078E-4 


6.68766E-8 








-2.33016E7 


4.24102E7 








-4.30546E7 


7.84699E7 









*Rows with the symbol * uses log A;, instead of Af, as y variables. 

**Rows with the symbol ** use rrienv = Msnv/M instead of -R as a; variables and 1/A as 2/ variables. 



Table 2. Fitting Results of A for Pop. II Stars 



Stage Mass(M0) A a 61 62 63 64 65 



1 


Aft 


1.06031 


0.08173 


-0.00436 


6.7847E-5 


-3.3429E-7 







A9 


0.57192 


0.03937 


-0.00221 


3.49192E-5 


-1.73046E-7 





2 


A6 


2.56108 


-0.75562 


0.1027 


-0.00495 


8.05436E-5 









1.41896 


-0.4266 


0.05792 


-0.00281 


4.61E-5 





3 


A6 


1.7814 


-0.17138 


0.00754 


-9.02652E-5 










A9 


0.99218 


-0.10082 


0.00451 


-5.53632E-5 








4 


At 


1.65914 


-0.10398 


0.0029 


-2.24862E-5 










At 


0.92172 


-0.06187 


0.00177 


-1.42677E-5 








5 


A6 


1.58701 


-0.06897 


0.00129 


-6.99399E-6 










Ag 


0.87647 


-0.04103 


7.91444E-4 


-4.41644E-6 








6 


A6 


1.527 


-0.04738 


6.1373E-4 


-2.36835E-6 










A9 


0.83636 


-0.02806 


3.73346E-4 


-1.47016E-6 








7 


At 


1.49995 


-0.03921 


4.2327E-4 


-1.37646E-6 










A9 


0.81688 


-0.02324 


2.5804E-4 


-8.54696E-7 








D 


A6 




-U.Uoio4 


Z.oOOZZIIj-4: 


7 fll 00017 7 


U 


u 




Ag 


0.79396 


-0.01903 


1.77574E-4 


-5.04262E-7 








9 


A6 


1.49196 


-0.03247 


3.08066E-4 


-9.53247E-7 










As 


0.805 


-0.02 


2.01872E-4 


-6.4295E-7 








10 


Af, 


1.73966 


-0.07074 


0.00178 


-1.73072E-5 










A9 


1.027 


-0.05685 


0.00164 


-1.65206E-5 








12 


A6 


1.63634 


-0.04646 


7.49351E-4 


-5.23622E-6 










Ag 


1.17934 


-0.08481 


0.00329 


-4.69096E-5 








14 


A6 


1.45573 


-0.00937 


-0.00131 


3.07004E-5 










Ag 


1.19526 


-0.08503 


0.00324 


-4.58919E-5 








16 


A6 


1.33378 


0.01274 


-0.00234 


4.6036E-5 










As 


1.17731 


-0.07834 


0.00275 


-3.58108E-5 








20 


A6 


1.27138 


0.00538 


-0.0012 


1.80776E-5 










Ag 


1.07496 


-0.05737 


0.00153 


-1.49005E-5 








1 


A6 


0.37294 


-0.05825 


0.00375 


-7.59191E-5 










Ag 


0.24816 


-0.04102 


0.0028 


-6.20419E-5 









Table 2 — Continued 



Stage Mass(M0) A a bi 62 63 64 65 



2 


2 




-103.92538 


25.37325 








-56.03478 


13.6749 


2 


3 




-12.40832 


1.59021 






As 


-6.47476 


0.8328 


2 


4 


A6 


-5.89253 


0.54296 






As 


-3.21299 


0.29583 


2 


5 


Ai, 


-0.67176 


0.07708 






As 


-0.38561 


0.0427 


2 


6 


Ab 


0.30941 


0.00965 






As 


0.14576 


0.00562 


2 


7 


A6 


0.44862 


0.00234 






As 


0.21873 


0.00154 


2 


8 




n 50221 


1 Qn21 F,-4 






As 


0.24748 


-9.9338E-5 


2 


9 


Ai, 


0.39342 


0.00259 






As 


0.20796 


6.62921E-4 


2 


10 


A6 


0.75746 


-0.00852 






As 


0.35355 


-0.00388 


2 


12 


Ai, 


0.09472 


1.59771 






As 


0.04621 


0.23145 


2 


14 


Ab 


0.08576 


0.50546 






As 


0.04015 


0.69617 


2 


16 


A6 


0.0774 


1.86228 






As 


0.03457 


0.73627 


2 


20 


Ai, 


0.05174 


1.52248 






As 


0.02187 


0.16728 


3 


1 


Ai, 


0.24012 


-0.01907 






As 


0.15504 


-0.01238 


3 


2 


A6 


0.5452 


0.00212 






As 


0.30594 


-9.58858E-4 



-2.03273 


0.0543 








-1.09533 


0.02925 








-0.06494 


8.69587E-4 








-0.03412 


4.58399E-4 








-0.01527 


1.38354E-4 








-0.00833 


7.55646E-5 








-0.00175 


1.1991E-5 








-9.6948E-4 


6.64455E-6 








■2.31975E-4 


1.26273E-6 








■1.30273E-4 


7.06459E-7 








■9.23152E-5 


4.67797E-7 








5.18806E-5 


2.60283E-7 










1 S072fiR-7 








■1.99272E-5 


9.47504E-8 








■4.97778E-5 


1.69533E-7 








■1.84663E-5 


6.58983E-8 








3.51646E-5 


-4.57725E-8 








1.56573E-5 


-1.98173E-8 








19.38009 


0.4662 


83.08091 





76.88504 


0.64758 


18.45577 





90.52062 


1.73793 


17.7433 





17.24357 


0.23011 


86.67988 





17.55952 


0.52273 


99.08578 





17.10365 


0.22597 


96.30501 





25.93257 


0.4196 


149.18541 





146.24645 


0.55973 








6.09529E-4 


-8.17819E-6 


4.83789E-8 


-1.04568E-10 


3.96633E-4 


-5.3329E-6 


3.16052E-8 


-6.84288E-11 


6.42941E-5 


-1.46783E-7 








1.12174E-4 


-1.04079E-6 


3.4564E-9 


-3.91536E-12 
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Table 2 — Continued 



Stage 


Mass(M0) 


A 


a 


6i 


62 


63 


&4 


65 


3 


3* 


Ab 


-0.475 


-0.00328 


1.31101E-4 


-6.03669E-7 


8.49549E-10 









As 


0.05434 


0.0039 


9.44609E-6 


-3.87278E-8 








3 


4* 


Ai, 


-0.2106 


-0.01574 


2.01107E-4 


-6.90334E-7 


7.92713E-10 









As 


0.36779 


-0.00991 


1.19411E-4 


-3.59574E-7 


3.33957E-10 





3 


5 


Ab 


-0.12027 


0.01981 


-2.27908E-4 


7.55556E-7 












As 


0.31252 


-0.00527 


3.60348E-5 


-3.22445E-8 








3 


6 


Ai, 


0.26578 


0.00494 


-7.02203E-5 


2.25289E-7 












As 


0.26802 


-0.00248 


6.45229E-6 


1.69609E-8 








3 


7 


A6 


0.8158 


-0.01633 


1.46552E-4 


-5.75308E-7 


8.77711E-10 









As 


0.26883 


-0.00219 


4.12941E-6 


1.33138E-8 








3 


8 


A(, 


0.74924 


-0.01233 


9.55715E-5 


-3.37117E-7 


4.67367E-10 









As 


0.25249 


-0.00161 


8.35478E-7 


1.25999E-8 








3 


9 


A(, 


0.73147 


-0.01076 


7.54308E-5 


-2.4114E-7 


2.95543E-10 









As 


0.31951 


-0.00392 


2.31815E-5 


-6.59418E-8 


7.99575E-11 





3 


10 


Ab 


-9.26519 


0.08064 


-2.30952E-4 


2.21986E-7 












As 


0.81491 


-0.00161 


-8.13352E-6 


1.95775E-8 








3 


12 


A(, 


-51.15252 


0.30238 


-5.95397E-4 


3.91798E-7 












As 


-13.44 


0.08141 


-1.641E-4 


1.106E-7 








3 


14 


A(, 


-140 


0.7126 


0.00121 


6.846E-7 












As 


-44.1964 


0.22592 


-3.85124E-4 


2.19324E-7 








3 


16 


Ab 


-358.4 


1.599 


-0.00238 


1.178E-6 












As 


-118.13757 


0.52737 


-7.8479E-4 


3.89585E-7 








3 


20 


Ab 


-436.00777 


1.41375 


-0.00153 


5.47573E-7 












As 


-144.53456 


0.46579 


-4.99197E-4 


1.78027E-7 









*Rows with the symbol * use logA^ instead of as y variables. 



